/****Replication DO-FILE for "Search for Spcies and Souls: Catholic Missions as Colonial State in the Philippines" by Dean Dulay, Singapore Management University*/

/**Run both the do-file and the data file on Desktop. An alternative is to manually**
load the data abd begin in line 8**/

use "replication_data_church.dta"

/**Main Outcomes: State Capacity and Development**/

reg govt_office catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo GO_OLS_1


reg govt_office catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo GO_OLS_2


reg govt_office catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times"
estadd local State  "$\times"
eststo GO_OLS_3


reg direct_tax_share_2010 catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo DTS_OLS_1

reg direct_tax_share_2010 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo DTS_OLS_2

reg direct_tax_share_2010 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo DTS_OLS_3

reg total_road_density catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo RD_OLS_1

reg total_road_density catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo RD_OLS_2

reg total_road_density catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo RD_OLS_3

reg schools_per_student_2015 catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo SPS_OLS_1

reg schools_per_student_2015 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo SPS_OLS_2

reg schools_per_student_2015 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo SPS_OLS_3


reg log_nightlights_2010 catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo NL_OLS_1

reg log_nightlights_2010 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo NL_OLS_2

reg log_nightlights_2010 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo NL_OLS_3

esttab GO_OLS_1 GO_OLS_2 GO_OLS_3 DTS_OLS_1 DTS_OLS_2 DTS_OLS_3 using "mainresults.tex", ///
prehead("\begin{tabular}{1*{6}{c}}\hline \hline") ///
fragment ///
mgroups("\multicolumn{3}{c}{Number of Government Offices} & \multicolumn{3}{c}{Direct Tax Shares}", span) title("Catholic Mission, State Capacity, and Development") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace
    
esttab RD_OLS_1 RD_OLS_2 RD_OLS_3 SPS_OLS_1 SPS_OLS_2 SPS_OLS_3 using "mainresults.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Road Density} & \multicolumn{3}{c}{Schools per Student}", span) title("Catholic Mission, State Capacity, and Development") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 

esttab NL_OLS_1 NL_OLS_2 NL_OLS_3 using "mainresults.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Night Light Luminosity}", span) title("Catholic Mission, State Capacity, and Development") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 


/************Mechanisms 1 and 2: Education and Urbanization****************/

reg schools_1918 catholic_mission pre_colonial sultanate i.provincefe if Population_1918>1000, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  ""
estadd local State ""
eststo EDUC_1

reg schools_1918 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if Population_1918>1000, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic "$\times$"
estadd local State ""
eststo EDUC_2

reg schools_1918 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if Population_1918>1000, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  "$\times$"
estadd local State "$\times$"
eststo EDUC_3

reg log_pop1903 catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  ""
estadd local State ""
eststo URB_1

reg log_pop1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic "$\times$"
estadd local State ""
eststo URB_2

reg log_pop1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  "$\times$"
estadd local State "$\times$"
eststo URB_3


esttab EDUC_1 EDUC_2 EDUC_3 URB_1 URB_2 URB_3 using "mechanisms1.tex", ///  
mgroups("\multicolumn{3}{c}{Schools in 1918} & \multicolumn{3}{c}{Population 1918}", span) title("Education and Urbanization Mechanisms") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace

/**Mechanisms 2: Occupational Specialization**/

reg agriculture_share_1903 catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo AS_OLS_1

reg agriculture_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo AS_OLS_2

reg agriculture_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo AS_OLS_3


reg manufacturing_share_1903 catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo MS_OLS_1

reg manufacturing_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo MS_OLS_2

reg manufacturing_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo MS_OLS_3


reg trade_share_1903 catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo TS_OLS_1

reg trade_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo TS_OLS_2

reg trade_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo TS_OLS_3

esttab AS_OLS_1 AS_OLS_2 AS_OLS_3 MS_OLS_1 MS_OLS_2 MS_OLS_3 TS_OLS_1 TS_OLS_2 TS_OLS_3 using "mechanisms2.tex", ///  
mgroups("\multicolumn{3}{c}{Agriculture Share} & \multicolumn{3}{c}{Manufacturing Share} & \multicolumn{3}{c}{Trade Share}", span) title("Structural Transformation") ///
drop (_cons pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 *.provincefe Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Secular State Controls" "N"))  ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace


/***Mechanism 3: Religious Competition***/

reg prot_mission catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo PM_OLS_1

reg prot_mission catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo PM_OLS_2

reg prot_mission catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo PM_OLS_3

reg prot_churches catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo PC_OLS_1

reg prot_churches catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo PC_OLS_2

reg prot_churches catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo PC_OLS_3

esttab PM_OLS_1 PM_OLS_2 PM_OLS_3 PC_OLS_1 PC_OLS_2 PC_OLS_3 using "mechanisms3.tex", ///  
mgroups("\multicolumn{3}{c}{Protestant Missions} & \multicolumn{3}{c}{Protestant Churches}", span) title("Religious Competition Mechanism") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Secular State Controls" "N"))  ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace


/*************ALTERNATIVE MECHANISMS**************/


reg telegraphlines_1903 catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo TL_OLS_1

reg telegraphlines_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo TL_OLS_2


reg telegraphlines_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo TL_OLS_3


esttab TL_OLS_1 TL_OLS_2 TL_OLS_3  using "alternative_mechanisms1.tex", ///  
mgroups("\multicolumn{3}{c}{Telegraph Lines}", span) title("Infrastructure") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Secular State Controls" "N"))  ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace


/**Variation Across Religious Orders**/


reg govt_office catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo GOJ_OLS_1


reg govt_office catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo GOJ_OLS_2


reg govt_office catholic_mission jesuit jesuit_mission_binary  pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times"
estadd local State  "$\times"
eststo GOJ_OLS_3


reg direct_tax_share_2010 catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo DTSJ_OLS_1

reg direct_tax_share_2010 catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo DTSJ_OLS_2

reg direct_tax_share_2010 catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo DTSJ_OLS_3

reg total_road_density catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo RDJ_OLS_1

reg total_road_density catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo RDJ_OLS_2

reg total_road_density catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo RDJ_OLS_3

reg schools_per_student_2015 catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo SPSJ_OLS_1

reg schools_per_student_2015 catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo SPSJ_OLS_2

reg schools_per_student_2015 catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo SPSJ_OLS_3


reg log_nightlights_2010 catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo NLJ_OLS_1

reg log_nightlights_2010 catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo NLJ_OLS_2

reg log_nightlights_2010 catholic_mission jesuit jesuit_mission_binary pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo NLJ_OLS_3

esttab GOJ_OLS_1 GOJ_OLS_2 GOJ_OLS_3 DTSJ_OLS_1 DTSJ_OLS_2 DTSJ_OLS_3 using "alternative_mechanisms2.tex", ///
prehead("\begin{tabular}{1*{6}{c}}\hline \hline") ///
fragment ///
mgroups("\multicolumn{3}{c}{Number of Government Offices} & \multicolumn{3}{c}{Direct Tax Shares}", span) title("Variation in Missionary Order") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace
    
esttab RDJ_OLS_1 RDJ_OLS_2 RDJ_OLS_3 SPSJ_OLS_1 SPSJ_OLS_2 SPSJ_OLS_3 using "alternative_mechanisms2.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Road Density} & \multicolumn{3}{c}{Schools per Student}", span) title("Variation in Missionary Order") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 

esttab NLJ_OLS_1 NLJ_OLS_2 NLJ_OLS_3 using "alternative_mechanisms2.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Night Light Luminosity}", span) title("Variation in Missionary Order") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 




/**Malaria**/

reg malaria_1903 catholic_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo ML_OLS_1


reg malaria_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo ML_OLS_2

reg malaria_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo ML_OLS_3


esttab ML_OLS_1 ML_OLS_2 ML_OLS_3 using "alternative_mechanisms3.tex", ///  
mgroups("\multicolumn{3}{c}{Malaria Deaths per 1,000}", span) title("Malaria Deaths per 1,000") ///
drop (_cons pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 *.provincefe Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Secular State Controls" "N"))  ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace


/**Friar Lands**/

reg govt_office mission_friar pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo GOFL_OLS_1


reg govt_office mission_friar pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo GOFL_OLS_2


reg govt_office mission_friar pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times"
estadd local State  "$\times"
eststo GOFL_OLS_3


reg direct_tax_share_2010 mission_friar pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo DTSFL_OLS_1

reg direct_tax_share_2010 mission_friar pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo DTSFL_OLS_2

reg direct_tax_share_2010 mission_friar pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo DTSFL_OLS_3

reg total_road_density mission_friar pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo RDFL_OLS_1

reg total_road_density mission_friar pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo RDFL_OLS_2

reg total_road_density mission_friar pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo RDFL_OLS_3

reg schools_per_student_2015 mission_friar pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm 
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo SPSFL_OLS_1

reg schools_per_student_2015 mission_friar pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo SPSFL_OLS_2

reg schools_per_student_2015 mission_friar pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo SPSFL_OLS_3


reg log_nightlights_2010 mission_friar pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo NLFL_OLS_1

reg log_nightlights_2010 mission_friar pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo NLFL_OLS_2

reg log_nightlights_2010 mission_friar pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo NLFL_OLS_3


esttab GOFL_OLS_1 GOFL_OLS_2 GOFL_OLS_3 DTSFL_OLS_1 DTSFL_OLS_2 DTSFL_OLS_3 using "alternative_mechanisms4.tex", ///
prehead("\begin{tabular}{1*{6}{c}}\hline \hline") ///
fragment ///
mgroups("\multicolumn{3}{c}{Number of Government Offices} & \multicolumn{3}{c}{Direct Tax Shares}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Secular State Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace
    
esttab RDFL_OLS_1 RDFL_OLS_2 RDFL_OLS_3 SPSFL_OLS_1 SPSFL_OLS_2 SPSFL_OLS_3 using "alternative_mechanisms4.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Road Density} & \multicolumn{3}{c}{Schools per Student}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Secular State Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 

esttab NLFL_OLS_1 NLFL_OLS_2 NLFL_OLS_3 using "alternative_mechanisms4.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Night Light Luminosity}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Secular State Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 



/**Summary Statistics****/

estpost summ catholic_mission govt_office direct_tax_share_2013 total_road_density schools_per_student_2015 log_nightlights_2010 pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda schools_1918 Population_1918 agriculture_share_1903 manufacturing_share_1903 trade_share_1903 prot_mission prot_churches telegraphlines_1903 malaria_1903 jesuit malaria catholic_mission_longevity
eststo summstats

esttab summstats using "Summary_Statistics.tex", cells("count(label(N)) mean(fmt(a3) label(Mean)) sd(fmt(a3) label(SD)) min(fmt(a3) label(Min)) max(fmt(a3) label(Max))") nomtitle nonumber noobs label wrap replace prehead(`"\begin{table}[htbp]\centering"' `"\setlength\tabcolsep{0pt}"' `"\caption{@title}"' `"\label{summary:stats}"' `"\begin{tabular}{lcccc}"' `"\toprule"') postfoot( ) title("Summary Statistics")




/***ROBUSTNESS***/

/**Mission Longevity**/

reg govt_office catholic_mission_longevity pre_colonial sultanate i.provincefe if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo GOL_OLS_1


reg govt_office catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo GOL_OLS_2


reg govt_office catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times"
estadd local State  "$\times"
eststo GOL_OLS_3


reg direct_tax_share_2010 catholic_mission_longevity pre_colonial sultanate i.provincefe if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Preclonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo DTSL_OLS_1

reg direct_tax_share_2010 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo DTSL_OLS_2

reg direct_tax_share_2010 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo DTSL_OLS_3

reg total_road_density catholic_mission_longevity pre_colonial sultanate i.provincefe if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo RDL_OLS_1

reg total_road_density catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo RDL_OLS_2

reg total_road_density catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo RDL_OLS_3

reg schools_per_student_2015 catholic_mission_longevity pre_colonial sultanate i.provincefe if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo SPSL_OLS_1

reg schools_per_student_2015 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Prcolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo SPSL_OLS_2

reg schools_per_student_2015 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.province if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo SPSL_OLS_3

reg log_nightlights_2010 catholic_mission_longevity pre_colonial sultanate i.provincefe if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo NLL_OLS_1

reg log_nightlights_2010 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo NLL_OLS_2

reg log_nightlights_2010 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if catholic_mission==1, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo NLL_OLS_3


esttab GOL_OLS_1 GOL_OLS_2 GOL_OLS_3 DTSL_OLS_1 DTSL_OLS_2 DTSL_OLS_3 using "catholic_longevity.tex", ///
prehead("\begin{tabular}{1*{6}{c}}\hline \hline") ///
fragment ///
mgroups("\multicolumn{3}{c}{Number of Government Offices} & \multicolumn{3}{c}{Direct Tax Shares}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Secular State Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace
    
esttab RDL_OLS_1 RDL_OLS_2 RDL_OLS_3 SPSL_OLS_1 SPSL_OLS_2 SPSL_OLS_3 using "catholic_longevity.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Road Density} & \multicolumn{3}{c}{Schools per Student}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 

esttab NLL_OLS_1 NLL_OLS_2 NLL_OLS_3 using "catholic_longevity.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Night Light Luminosity}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 


/**Main Results by Century**/

/*****1500s mission*******/

reg govt_office c1500_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Pre-Colonial  "$\times$"
estadd local Geographic  ""
estadd local Secular State  ""
eststo GO15_OLS_1


reg govt_office c1500_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo GO15_OLS_2


reg govt_office c1500_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times"
estadd local State  "$\times"
eststo GO15_OLS_3


reg direct_tax_share_2010 c1500_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo DTS15_OLS_1

reg direct_tax_share_2010 c1500_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo DTS15_OLS_2

reg direct_tax_share_2010 c1500_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo DTS15_OLS_3

reg total_road_density c1500_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo RD15_OLS_1

reg total_road_density c1500_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo RD15_OLS_2

reg total_road_density c1500_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo RD15_OLS_3

reg schools_per_student_2015 c1500_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo SPS15_OLS_1

reg schools_per_student_2015 c1500_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo SPS15_OLS_2

reg schools_per_student_2015 c1500_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo SPS15_OLS_3


reg log_nightlights_2010 c1500_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo NL15_OLS_1

reg log_nightlights_2010 c1500_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo NL15_OLS_2

reg log_nightlights_2010 c1500_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo NL15_OLS_3


esttab GO15_OLS_1 GO15_OLS_2 GO15_OLS_3 DTS15_OLS_1 DTS15_OLS_2 DTS15_OLS_3 using "catholic_15.tex", ///
prehead("\begin{tabular}{1*{6}{c}}\hline \hline") ///
fragment ///
mgroups("\multicolumn{3}{c}{Number of Government Offices} & \multicolumn{3}{c}{Direct Tax Shares}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace
    
esttab RD15_OLS_1 RD15_OLS_2 RD15_OLS_3 SPS15_OLS_1 SPS15_OLS_2 SPS15_OLS_3 using "catholic_15.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Road Density} & \multicolumn{3}{c}{Schools per Student}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 

esttab NL15_OLS_1 NL15_OLS_2 NL15_OLS_3 using "catholic_15.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Night Light Luminosity}", span) title("Catholic Missions in 1500") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 


/********16th Century******/

reg govt_office c1600_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo GO16_OLS_1


reg govt_office c1600_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo GO16_OLS_2


reg govt_office c1600_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times"
estadd local State  "$\times"
eststo GO16_OLS_3


reg direct_tax_share_2010 c1600_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo DTS16_OLS_1

reg direct_tax_share_2010 c1600_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo DTS16_OLS_2

reg direct_tax_share_2010 c1600_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo DTS16_OLS_3

reg total_road_density c1600_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo RD16_OLS_1

reg total_road_density c1600_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo RD16_OLS_2

reg total_road_density c1600_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo RD16_OLS_3

reg schools_per_student_2015 c1600_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo SPS16_OLS_1

reg schools_per_student_2015 c1600_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo SPS16_OLS_2

reg schools_per_student_2015 c1600_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo SPS16_OLS_3


reg log_nightlights_2010 c1600_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo NL16_OLS_1

reg log_nightlights_2010 c1600_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo NL16_OLS_2

reg log_nightlights_2010 c1600_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo NL16_OLS_3


esttab GO16_OLS_1 GO16_OLS_2 GO16_OLS_3 DTS16_OLS_1 DTS16_OLS_2 DTS16_OLS_3 using "catholic_16.tex", ///
prehead("\begin{tabular}{1*{6}{c}}\hline \hline") ///
fragment ///
mgroups("\multicolumn{3}{c}{Number of Government Offices} & \multicolumn{3}{c}{Direct Tax Shares}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace
    
esttab RD16_OLS_1 RD16_OLS_2 RD16_OLS_3 SPS16_OLS_1 SPS16_OLS_2 SPS16_OLS_3 using "catholic_16.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Road Density} & \multicolumn{3}{c}{Schools per Student}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 

esttab NL16_OLS_1 NL16_OLS_2 NL16_OLS_3 using "catholic_16.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Night Light Luminosity}", span) title("Catholic Missions in 1600") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 


/*******1700s Mission********/

reg govt_office c1700_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo GO17_OLS_1


reg govt_office c1700_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo GO17_OLS_2


reg govt_office c1700_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times"
estadd local State  "$\times"
eststo GO17_OLS_3


reg direct_tax_share_2010 c1700_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo DTS17_OLS_1

reg direct_tax_share_2010 c1700_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo DTS17_OLS_2

reg direct_tax_share_2010 c1700_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo DTS17_OLS_3

reg total_road_density c1700_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo RD17_OLS_1

reg total_road_density c1700_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo RD17_OLS_2

reg total_road_density c1700_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo RD17_OLS_3

reg schools_per_student_2015 c1700_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo SPS17_OLS_1

reg schools_per_student_2015 c1700_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo SPS17_OLS_2

reg schools_per_student_2015 c1700_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo SPS17_OLS_3

reg log_nightlights_2010 c1700_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo NL17_OLS_1

reg log_nightlights_2010 c1700_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo NL17_OLS_2

reg log_nightlights_2010 c1700_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo NL17_OLS_3


esttab GO17_OLS_1 GO17_OLS_2 GO17_OLS_3 DTS17_OLS_1 DTS17_OLS_2 DTS17_OLS_3 using "catholic_17.tex", ///
prehead("\begin{tabular}{1*{6}{c}}\hline \hline") ///
fragment ///
mgroups("\multicolumn{3}{c}{Number of Government Offices} & \multicolumn{3}{c}{Direct Tax Shares}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace
    
esttab RD17_OLS_1 RD17_OLS_2 RD17_OLS_3 SPS17_OLS_1 SPS17_OLS_2 SPS17_OLS_3 using "catholic_17.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Road Density} & \multicolumn{3}{c}{Schools per Student}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 

esttab NL17_OLS_1 NL17_OLS_2 NL17_OLS_3 using "catholic_17.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Night Light Luminosity}", span) title("Catholic Missions in 1600") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 


/*****1800 Mission*******/

reg govt_office c1800_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo GO18_OLS_1


reg govt_office c1800_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo GO18_OLS_2


reg govt_office c1800_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times"
estadd local State  "$\times"
eststo GO18_OLS_3


reg direct_tax_share_2010 c1800_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo DTS18_OLS_1

reg direct_tax_share_2010 c1800_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo DTS18_OLS_2

reg direct_tax_share_2010 c1800_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo DTS18_OLS_3

reg total_road_density c1800_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo RD18_OLS_1

reg total_road_density c1800_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo RD18_OLS_2

reg total_road_density c1800_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo RD18_OLS_3

reg schools_per_student_2015 c1800_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo SPS18_OLS_1

reg schools_per_student_2015 c1800_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo SPS18_OLS_2

reg schools_per_student_2015 c1800_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo SPS18_OLS_3


reg log_nightlights_2010 c1800_mission pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo NL18_OLS_1

reg log_nightlights_2010 c1800_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo NL18_OLS_2

reg log_nightlights_2010 c1800_mission pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo NL18_OLS_3

esttab GO18_OLS_1 GO18_OLS_2 GO18_OLS_3 DTS18_OLS_1 DTS18_OLS_2 DTS18_OLS_3 using "catholic_18.tex", ///
prehead("\begin{tabular}{1*{6}{c}}\hline \hline") ///
fragment ///
mgroups("\multicolumn{3}{c}{Number of Government Offices} & \multicolumn{3}{c}{Direct Tax Shares}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace
    
esttab RD18_OLS_1 RD18_OLS_2 RD18_OLS_3 SPS18_OLS_1 SPS18_OLS_2 SPS18_OLS_3 using "catholic_18.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Road Density} & \multicolumn{3}{c}{Schools per Student}", span) title("Friar Lands") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 

esttab NL18_OLS_1 NL18_OLS_2 NL18_OLS_3 using "catholic_18.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Night Light Luminosity}", span) title("Catholic Missions in 1600") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "SGovernment Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 


/****No Outliers****/


reg govt_office catholic_mission pre_colonial sultanate i.provincefe if govt_office!=0, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo GO_OLS_1_NO


reg govt_office catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if govt_office!=0, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo GO_OLS_2_NO


reg govt_office catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if govt_office!=0, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times"
estadd local State  "$\times"
eststo GO_OLS_3_NO


reg direct_tax_share_2010 catholic_mission pre_colonial sultanate i.provincefe if direct_tax_share_2010>.0023492 & direct_tax_share_2010<.2910103, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo DTS_OLS_1_NO

reg direct_tax_share_2010 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if direct_tax_share_2010>.0023492 & direct_tax_share_2010<.2910103, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo DTS_OLS_2_NO

reg direct_tax_share_2010 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if direct_tax_share_2010>.0023492 & direct_tax_share_2010<.2910103, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo DTS_OLS_3_NO

reg total_road_density catholic_mission pre_colonial sultanate i.provincefe if total_road_density>0 & total_road_density<3.140251, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo RD_OLS_1_NO

reg total_road_density catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if total_road_density>0 & total_road_density<3.140251, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo RD_OLS_2_NO

reg total_road_density catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if total_road_density>0 & total_road_density<3.140251, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo RD_OLS_3_NO

reg schools_per_student_2015 catholic_mission pre_colonial sultanate i.provincefe if schools_per_student_2015>.2038437 & schools_per_student_2015<1.415167, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo SPS_OLS_1_NO

reg schools_per_student_2015 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if schools_per_student_2015>.2038437 & schools_per_student_2015<1.415167, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo SPS_OLS_2_NO

reg schools_per_student_2015 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if schools_per_student_2015>.2038437 & schools_per_student_2015<1.415167, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo SPS_OLS_3_NO


reg log_nightlights_2010 catholic_mission pre_colonial sultanate i.provincefe if log_nightlights_2010>-2.429473 & log_nightlights_2010<3.005168, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo NL_OLS_1_NO

reg log_nightlights_2010 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if log_nightlights_2010>-2.429473 & log_nightlights_2010<3.005168, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo NL_OLS_2_NO

reg log_nightlights_2010 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if log_nightlights_2010>-2.429473 & log_nightlights_2010<3.005168, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo NL_OLS_3_NO



esttab GO_OLS_1_NO GO_OLS_2_NO GO_OLS_3_NO DTS_OLS_1_NO DTS_OLS_2_NO DTS_OLS_3_NO using "mainresults_nooutliers.tex", ///
prehead("\begin{tabular}{1*{6}{c}}\hline \hline") ///
fragment ///
mgroups("\multicolumn{3}{c}{Number of Government Offices} & \multicolumn{3}{c}{Direct Tax Shares}", span) title("Catholic Mission, State Capacity, and Development (No Outliers)") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace
    
esttab RD_OLS_1_NO RD_OLS_2_NO RD_OLS_3_NO SPS_OLS_1_NO SPS_OLS_2_NO SPS_OLS_3_NO using "mainresults_nooutliers.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Road Density} & \multicolumn{3}{c}{Schools per Student}", span) title("Catholic Mission, State Capacity, and Development (No Outliers)") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 

esttab NL_OLS_1_NO NL_OLS_2_NO NL_OLS_3_NO using "mainresults_nooutliers.tex", ///
fragment ///
append ///
mgroups("\multicolumn{3}{c}{Night Light Luminosity}", span) title("Catholic Mission, State Capacity, and Development (No Outliers)") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N")) ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se 


/****************Appendix for Mechanisms******************/


/****Education and Urbanization (Longevity)******/


reg schools_1918 catholic_mission_longevity pre_colonial sultanate i.provincefe if Population_1918>1000, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  ""
estadd local State ""
eststo EDUC_1_LONG

reg schools_1918 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if Population_1918>1000, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic "$\times$"
estadd local State ""
eststo EDUC_2_LONG

reg schools_1918 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if Population_1918>1000, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  "$\times$"
estadd local State "$\times$"
eststo EDUC_3_LONG

reg log_pop1903 catholic_mission_longevity pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  ""
estadd local State ""
eststo URB_1_LONG

reg log_pop1903 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic "$\times$"
estadd local State ""
eststo URB_2_LONG

reg log_pop1903 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  "$\times$"
estadd local State "$\times$"
eststo URB_3_LONG


esttab EDUC_1_LONG EDUC_2_LONG EDUC_3_LONG URB_1_LONG URB_2_LONG URB_3_LONG using "mechanisms1_longevity.tex", ///  
mgroups("\multicolumn{3}{c}{Schools in 1918} & \multicolumn{3}{c}{Population 1918}", span) title("Education and Urbanization Mechanisms (Longevity Measure)") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N"))  ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace



/***Education and Urbanization (No Outliers)***/

reg schools_1918 catholic_mission pre_colonial sultanate i.provincefe if Population_1918>1000 & schools_1918>1 & schools_1918<13, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  ""
estadd local State ""
eststo EDUC_1_NO

reg schools_1918 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if Population_1918>1000 & schools_1918>1 & schools_1918<13, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic "$\times$"
estadd local State ""
eststo EDUC_2_NO

reg schools_1918 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if Population_1918>1000 & schools_1918>1 & schools_1918<13, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  "$\times$"
estadd local State "$\times$"
eststo EDUC_3_NO

reg log_pop1903 catholic_mission pre_colonial sultanate i.provincefe if log_pop1903>7.17549 & log_pop1903<9.80444, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  ""
estadd local State ""
eststo URB_1_NO

reg log_pop1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if log_pop1903>7.17549 & log_pop1903<9.80444, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic "$\times$"
estadd local State ""
eststo URB_2_NO

reg log_pop1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if log_pop1903>7.17549 & log_pop1903<9.80444, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial "$\times$"
estadd local Geographic  "$\times$"
estadd local State "$\times$"
eststo URB_3_NO


esttab EDUC_1_NO EDUC_2_NO EDUC_3_NO URB_1_NO URB_2_NO URB_3_NO using "mechanisms1_nooutliers.tex", ///  
mgroups("\multicolumn{3}{c}{Schools in 1918} & \multicolumn{3}{c}{Population 1918}", span) title("Education and Urbanization Mechanisms (No Outliers)") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N"))  ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace




/*****Structural Transformation (Longevity)****/

reg agriculture_share_1903 catholic_mission_longevity pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo AS_OLS_1_LONG

reg agriculture_share_1903 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo AS_OLS_2_LONG

reg agriculture_share_1903 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo AS_OLS_3_LONG


reg manufacturing_share_1903 catholic_mission_longevity pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo MS_OLS_1_LONG

reg manufacturing_share_1903 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo MS_OLS_2_LONG

reg manufacturing_share_1903 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo MS_OLS_3_LONG


reg trade_share_1903 catholic_mission_longevity pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo TS_OLS_1_LONG

reg trade_share_1903 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo TS_OLS_2_LONG

reg trade_share_1903 catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo TS_OLS_3_LONG

esttab AS_OLS_1_LONG AS_OLS_2_LONG AS_OLS_3_LONG MS_OLS_1_LONG MS_OLS_2_LONG MS_OLS_3_LONG TS_OLS_1_LONG TS_OLS_2_LONG TS_OLS_3_LONG using "mechanisms2_longevity.tex", ///  
mgroups("\multicolumn{3}{c}{Agriculture Share} & \multicolumn{3}{c}{Manufacturing Share} & \multicolumn{3}{c}{Trade Share}", span) title("Structural Transformation (Longevity)") ///
drop (_cons pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 *.provincefe Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N"))  ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace


/****Occupational Specialization (No Outliers)**/


reg agriculture_share_1903 catholic_mission pre_colonial sultanate i.provincefe if agriculture_share_1903>.0409792 & agriculture_share_1903<.3300832, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo AS_OLS_1_NO

reg agriculture_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if agriculture_share_1903>.0409792 & agriculture_share_1903<.3300832, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo AS_OLS_2_NO

reg agriculture_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if agriculture_share_1903>.0409792 & agriculture_share_1903<.3300832, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo AS_OLS_3_NO


reg manufacturing_share_1903 catholic_mission pre_colonial sultanate i.provincefe if manufacturing_share_1903>.0124732 & manufacturing_share_1903<.2919544, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo MS_OLS_1_NO

reg manufacturing_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if manufacturing_share_1903>.0124732 & manufacturing_share_1903<.2919544, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo MS_OLS_2_NO

reg manufacturing_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if manufacturing_share_1903>.0124732 & manufacturing_share_1903<.2919544, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo MS_OLS_3_NO


reg trade_share_1903 catholic_mission pre_colonial sultanate i.provincefe if trade_share_1903>.0018018 & trade_share_1903<.0844196, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo TS_OLS_1_NO

reg trade_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp if trade_share_1903>.0018018 & trade_share_1903<.0844196, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo TS_OLS_2_NO

reg trade_share_1903 catholic_mission pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda if trade_share_1903>.0018018 & trade_share_1903<.0844196, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times$"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo TS_OLS_3_NO

esttab AS_OLS_1_NO AS_OLS_2_NO AS_OLS_3_NO MS_OLS_1_NO MS_OLS_2_NO MS_OLS_3_NO TS_OLS_1_NO TS_OLS_2_NO TS_OLS_3_NO using "mechanisms2_nooutliers.tex", ///  
mgroups("\multicolumn{3}{c}{Agriculture Share} & \multicolumn{3}{c}{Manufacturing Share} & \multicolumn{3}{c}{Trade Share}", span) title("Structural Transformation (No Outliers)") ///
drop (_cons pre_colonial sultanate SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 *.provincefe Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N"))  ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace



/*****Religious Competition (Longevity)******/


reg prot_mission catholic_mission_longevity pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo PM_OLS_1_LONG

reg prot_mission catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo PM_OLS_2_LONG

reg prot_mission catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo PM_OLS_3_LONG

reg prot_churches catholic_mission_longevity pre_colonial sultanate i.provincefe, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  ""
estadd local State  ""
eststo PC_OLS_1_LONG

reg prot_churches catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  ""
eststo PC_OLS_2_LONG

reg prot_churches catholic_mission_longevity pre_colonial sultanate i.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda, cluster(PROVINCE) robust
estadd ysumm
estadd local FE "$\times"
estadd local Precolonial  "$\times$"
estadd local Geographic  "$\times$"
estadd local State  "$\times$"
eststo PC_OLS_3_LONG

esttab PM_OLS_1_LONG PM_OLS_2_LONG PM_OLS_3_LONG PC_OLS_1_LONG PC_OLS_2_LONG PC_OLS_3_LONG using "mechanisms3_longevity.tex", ///  
mgroups("\multicolumn{3}{c}{Protestant Missions} & \multicolumn{3}{c}{Protestant Churches}", span) title("Religious Competition Mechanism") ///
drop (_cons pre_colonial sultanate *.provincefe SQ1 SQ2 SQ3 SQ4 SQ5 SQ7 Coast_miles TRI_miles Precip Temp Manila_miles encomienda) stat(ymean FE Precolonial Geographic State N, labels("Mean of DV" "Province FE" "Pre-Colonial Controls" "Geographic Controls" "Government Controls" "N"))  ///
label nomtitles star(* 0.10 ** 0.05 *** 0.01) se replace





























